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The results of a theoretical investigation of the energy loss of charged particles in a magnetized classical 
plasma due to the electric field fluctuations are reported. The energy loss for a test particle is calculated through 
the linear-response theory. At vanishing magnetic field the electric field fluctuations lead to an energy gain of 
the charged particle for all velocities. It has been shown that in the presence of strong magnetic field this effect 
' occurs only at low-velocities. In the opposite case of high-velocities the test particle systematically loses its 

, energy due to the interaction with a stochastic electric field. The net effect of the fluctuations is the systematic 

reduction of the total energy loss (i.e. the sum of the polarization and stochastic energy losses) at vanishing 
magnetic field and reduction or enhancement at strong field depending on the velocity of the particle. It is found 
that the energy loss of the slow heavy ion contains an anomalous term which depends logarithmically on the 
projectile mass. The physical origin of this anomalous term is the coupling between the cyclotron motion of 
the plasma electrons and the long-wavelength, low-frequency fluctuations produced by the projectile ion. This 
effect may strongly enhance the stochastic energy gain of the particle. 
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I. INTRODUCTION 



C/2 

Q \ There is an ongoing interest in the theory of interaction of charged particle beams with plasmas. Although most theoretical 
" ^ ■ works have reported on the energy loss of ions in a plasma without magnetic field, the strongly magnetized case has not yet 
received as much attention as the field-free case. The energy loss of ion beams and the related processes in magnetized plasmas 
(~| are important in many areas of physics such as transport, heating, magnetic confinement of thermonuclear plasmas, and astro- 
CL( physics. Recent applications are the cooling of heavy-ion beams by electrons OJ-Ill] and the energy transfer for magnetized target 
fusion [4], or heavy-ion inertial confinement fusion. 

For a theoretical description of the energy loss of ions in a plasma, there exist some standard approaches. The dielectric linear 
^ response (LR) treatment considers the ion as a perturbation of the target plasma and the stopping is caused by the polarization 
00 of the surrounding medium. It is generally valid if the ion couples weakly to the target. Since the early 1960s, a number of 
calculations of the stopping power within LR treatment in a magnetized plasma have been presented (see Refs. ll5l-[T6ll and 
OO , references therein). Alternatively, the stopping is calculated as a result of the energy transfers in successive binary collisions 
' (BCs) between the ion and the electrons LI 7,-20.1 . Here it is necessary to consider appropriate approximations for the screening 
. of the Coulomb potential by the plasma H- However, significant gaps between these approaches involve the crucial ion stopping 
along magnetic field B and perpendicular to it. In particular, at high B values, the EC predicts a vanishingly parallel energy loss, 
which remains at variance with the nonzero LR one. Also challenging BC-LR discrepancies persist in the transverse direction, 
[ especially for vanishingly small ion projectile velocity v when the friction coefficient contains an anomalous term diverging 
J> logarithmically at v — > iflolfTTll . For calculation of the energy loss of an ion two new alternative approaches have been recently 
suggested. One of these methods is specifically aimed at a low-velocity energy loss, which is expressed in terms of velocity- 
velocity correlation and, hence, to a diffusion coefficient |16]. Next in Ref. [12] using the Bhatnagar-Krook-Gross approach 
based on the Boltzmann-Poisson equations for a collisional and magnetized classical plasma the energy loss of an ion is studied 
through a LR approach, which is constructed such that it conserves particle local number. 

Previous treatments of the energy loss of charged particles in a magnetized plasma J^-IMl were obtained by treating the 
plasma in an average manner, i.e., microscopic fluctuations were neglected. However, any media and, in particular, any plasma 
being a statistical system, is characterized by stochastic fluctuations. These fluctuations couple to external perturbations and the 
response of the medium to these perturbations can be expressed through suitable correlation functions of the microscopically 
fluctuating variables. It is well known that the motion of charged particles in such an environment is stochastic in nature and 
resembles Brownian motion. The effect of the electromagnetic field fluctuations on the energy loss of a charged particle moving 



•Electronic address: 
tElectronic address: 



hrachya@irphe.am | 



claude.deutsch@pgp.u-psud.fr] 



2 



in an unmagnetized classical plasma has been studied by several authors in the past ||2T1 - |28| . This effect leads to an energy gain 
of the particles and is most effective in the low-velocity limit. In addition the effect of the fluctuations are mainly important 
for light particles (e.g., for an electron) since this effect is inversely proportional to the mass of the particle and is negligible for 
heavy ions. Recently similar problem has been formulated for the energy loss in a so-called quark-gluon plasma [29] which is 
expected to be formed in relativistic heavy-ion collision experiments. Since the subject of the energy loss is of topical interest, it 
is the principal motivation of the present paper to quantitatively investigate the effect of microscopic fluctuations on the energy 
loss of a charged particle passing through an equilibrium, weakly-coupled magnetized plasma. In particular, our main objective 
is to study through the LR theory the effect of a strong magnetic field on the energy loss process arising due to the fluctuations 
in a plasma. The sequel is structured as follows. Interactions of a gyrating projectile with target stochastic electric fields are 
investigated at length in Sec. Resulting energy loss rates are respectively explained in Sec. |III]for a target plasma with no 
magnetic field and then with an infinite one as well as in Sec.|IV]for a finite magnetic field and at low-velocity of the particle. 



II. INTERACTION OF A GYRATING PARTICLE WITH PLASMA STOCHASTIC ELECTRIC FIELD 

The dielectric function approach was adopted earlier to calculate the energy loss of charged particles due to polarization effects 
of the medium It is assumed that the energy lost by the particle is small compared to the energy of the particle itself 

so that the change in the velocity of the particle during the motion may be neglected, i.e, the particle moves in a straight-line 
(in S = case) or helical (in S ^ case) trajectory. We consider a nonrelativistic projectile particle with charge Ze (— e is 
the electron charge) and with a velocity v(i), which moves in a magnetized classical plasma with a constant magnetic field B. 
We ignore any role of the electron spin or magnetic moment due to the nonrelativistic motion of the particle and the plasma 
electrons. The energy loss of a particle is determined by the work of the retarding forces acting on the particle in the plasma 
from the electric field generate d by the particle itself while moving. So the energy loss of the particle per unit time, i.e. energy 
loss rate (ELR), is given by ll2l-tl6ll 

^ = Zev(i).E(r(t),t), (1) 



where the electric field is taken at the location r(t) of the test particle. Within linear response theory the correlation function 
of the fluctuations of charge and current densities and the electromagnetic fields in the medium with space-time dispersion are 
completely determined in terms of the dielectric tensor of the medium. The total electric field E(x, t) induced in the plasma can 
be related to the external charge of the test particle by solving Poisson's equation (see, e.g., Ref. ll28ll ) 

{2nf J-00 J k^e{k,Lo)J_^ 

Here e(k,aj) — {kikj /k'^)eij{'k,uj) is the longitudinal dielectric function and ey (k,a;) is the dielectric tensor of the plasma, 
respectively. 

The previous formula for the energy loss in Eq. ([T) does not take into account the fluctuations of the electric fields in the 
plasma and the particle recoil in interactions. To accommodate these effects it is necessary to replace Eq. ([TJ with L24i,,25i1 . 

^^Ze{^{t)-£{rit),t)), (3) 

where (• • • ) denotes the statistical averaging operation. It is to be noted that here two averaging procedures are performed: (i) 
an ensemble average with respect to the equilibrium distribution function and (ii) a time average over random fluctuations in 
plasma. These two operations are commuting and only after both of them are performed the average quantity takes up a smooth 
value |26]. In the following, we will explicitly denote the ensemble average by (• • • )a wherever required to avoid possible 
confusion. 

The electric field £(x, t) = E(x, t) + E(x, t) in Eq. (O consists of the induced field E(x, t) given by Eq. ^ and a sponta- 
neously generated microscopic field E(x, t), the latter being a random function of position and time. 

The classical equation of motion of the particle with charge Ze and mass M moving in the fluctuating electric field 5(x, t) 
and a uniform magnetic field B have the form 

v(t) = [v (t) xh] + ^£ {r{t),t) , (4) 



where — \Z\/Z, ilc = \Z\eB/Mc is the cyclotron frequency of the test particle, b = B/B is the unit vector along the 
(external) magnetic field, £{x, t) is the total electric field involved in Eq. (|3]l, r(t) is the coordinate of the particle. Note that 
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Eq. (|4|i is the generalization of the Langevin equation in the presence of a magnetic field where the role of the friction force plays 
the polarization force determined by a regular electric field E(x, t) given by Eq. (|2]). 

Introducing the cyclotron rotation matrix Cij (t) the equation of motion (|4]i has a formal solution 

v^ (t) = V - (t) + — / dt'TZ,, {t - t') £, {r{t'),t') , (5) 

J —OO 

X, (t) = X, (t) + J dt'C,, {t - t') E, (r(i'), t') , (6) 

where TZij{t) = {1 /D,c)dCij{t) / dt, R(i) and V(i) with the components Xi{t) and Viit), respectively, are the initial (at t = 
—oo) unperturbed coordinate and velocity of the particle free cyclotron motion, respectively 

Xi (<) = XQi + -^^tj {t) voj, (7) 

v,{t)^n^j{t)voj. 

Here the variables vq and Tq with the components voi and xoi, respectively, are independent and are determined by initial 
conditions for the free cyclotron motion. Also 

Cij (t) = hihjVtct + {Sij - bibj) sin (ficO + <;eijsbs [1 - cos (flat)] , (8) 

where hi is the component of the unit vector b and e^^ is the fully antisymmetrical unit tensor of rank three. 

The mean change of the energy of the particle per unit time is given by Eq. (|3]l. The stochastic time dependence as embodied in 
Eq. (|3]l comes because of explicit fluctuations in time and motion of the particle from one field point to another The dependence 
on the latter will be expanded about the free cyclotron motion where the former will be left untouched. This is done by taking 
a time interval At sufficiently large with respect to the time scale of random electric field fluctuations in the plasma (T2) but 
small compared with the time during which the particle motion changes appreciably (ri), T2 <C At <C ti. Explicitly, ri ^ 
[e2r3/2ln(l/e)]"i and T2 (eT^/^^-i |[2^, where e = {AnneX^)^^ is the plasma parameter (Ad = {T/ATrriee^Y^^ and Ue 
are the Debye screening length and electron density, respectively) and T is the plasma temperature measured in energy units. Of 
course, in the present context of a magnetized plasma the time scales ti and T2 may be strongly modified by the magnetic field, 
see, e.g., Ref. 11211. In fact in the Coulomb logarithm determining the time ti, the Debye length is replaced by the cyclotron 
radius [3(J]. This results in a logarithmic dependence on magnetic field in the relaxation time ti. On the other hand, at strong 
magnetic fields T2 is replaced by the cyclotron period T2 ^ <^c^^, where ujc = eB/mc is the cyclotron frequency of the electrons. 
Thus, independently of the strength of the magnetic field and in the weak coupling limit with e <C 1, these scales are widely 
separated. In Eqs. (|5]l and (|6]l keeping only the leading order terms in the expansion we obtain 

(t) =V,{t) + — J dt'n^, {t - t') £, [t') , (9) 



(r{t),t) = it) + f dt'C.s {t ~ t')£, [t') , (10) 

where £{t) = 5(R(t), t) and R(t) is determined by Eq. (|7]l. Substituting Eqs. (|9]l and (fTOl i into Eq. ((Sjl and neglecting the third 
order correlations we obtain 

dt dt dt ' 

with 

dWQ 



dt 



- Ze\{t) ■ {£ {t)) , (12) 



Because the mean value of the fluctuating part of the electric field equals zero, (E(x, t)) =0, {£ (t)) equals the electric field 
E(R(t) , t) produced by the particle itself in the plasma. Therefore the first term in Eq. (fill determined by Eq. (ITZt corresponds 
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to the usual polarization energy loss of the particle calculated for instance in Ref. ||9[| fo r a magnetized plasma. For a gyrating 
projectile particle Eq. (fTZt has been evaluated explicitly in Ref. |9] (see also Ref. 111511 ). Averaging the energy loss dWo/dt 
with respect to the test particle cyclotron period 27r/17c (we denote this quantity as {dWo / dt) = Sq) we obtain (see Ref. |@] for 



details) 

72„2 °o 



^0 = ^j:J ^'^^^"(^-)i-7(k:c::M' '''^ 



where (nO^) — k\\Vo\\ + nilc, fc|| = (k • b), k'^ — k'^ — fc^, J„ is the Bessel function of the nth order Here vo\\ and vq± ^ are 
the unperturbed test particle velocities parallel and transverse to the magnetic field direction b, respectively, and a — uo±/^^c is 
the cyclotron radius. 

Assuming that there exists a hierarchy of scales T2 At ti, it can be shown that the polarization field does not contribute 
to leading order in the correlations functions appearing in the second term in Eq. (fTTT i as in the case of higher-order Fokker- 
Planck coefficients i22[l28l] . These terms correspond to the statistical change in the energy of the projectile particle in the plasma 
due to the fluctuations of the electric field as well as the velocity of the particle under the influence of this field. The first term 
in Eq. (flJl i corresponds to the statistical part of the dynamic friction due to the space-time correlation in the fluctuations in 
the electric field, whereas the second one corresponds to the average change in the energy of the moving particle due to the 
correlation between the fluctuation in the velocity of the particle and the fluctuation in the electrical field in the plasma. The 
temporal averaging in Eq. (fT3] l, by definition, includes many random fluctuations over the mean cyclotron motion. However, 
the correlation function of these fluctuations are suppressed beyond their characteristic time scales. This allows us to formally 
extend the limits in time integrations in Eq. ST3[ to infinity. 

Consider now the energy loss rate ( fTJt of the particle due to the interaction with stochastic plasma electric field. We consider 
homogeneous plasma and only electrostatic electric fields. In this case the correlation function (£i(x, t)£j{x' , t'))^ in Eq. ( fT3] l is 
a function of x — x' and t — t' . Introducing the power spectrum /'(k, uj) (which is determined by (x, t))a) of the fluctuating 
electric field this correlation function reads ll24ll25ll 

(£,(x,0f,(x',O)a- / P ^P(k,..). (15) 

Similar expression can be obtained for the correlation function involved in the second term of Eq. ( fTST l and containing the field 
gradient. Substituting these expressions into Eq. (fT?t we arrive at 



dt MVlc .1 k^ ./-n., ' ' V dt 

where 



dujP(k,uj)[—+iu]H{k,Lo,t), (16) 



H{k,uj,t) = f dt'C{t^t')e'^<'^^'^-'^(^'">'>"^^'~''\ 

— OO 



(17) 



C (t) = hk-jCij {t) = k^Vtct + k\ sin {Vtct) . (18) 

For evaluation of the function _ff(k, we substitute the unperturbed coordinate (see Eq. (|7)) of the free helical motion 
of the particle into Eq. ( [TtI i. The time-integration in Eq. ([TtI i can be performed using the Fourier series representation of the 
exponential function 1.31.] . After straightforward integration we obtain 

oo ( —k^il 

H{k,io,t)= V J„(fcia)e^("-')^<^ ^ -^Jiik^a) (19) 

[[u-Ci (k) - «0] 



n,l — ~oo 

kl 1 

'~2uj~Ci (k) - iO 



[Ji+i {k^a)e-'^ -Ji-i{k^a)e'^] 



Here V' — ^ ^ ^ct, ip is the initial phase of the cyclotron motion of the test particle with vq^ — vo± coscp and 
'^Oy = "^O-L sin if, 9 is the azimuthal angle of the transverse wave vector kj^ with k^ — k± cos 9 and ky — k± sin 9. Next 
we substitute Eq. ( fT9] l into Eq. (fT6] l and average the energy loss rate dWst / dt with respect to the particle cyclotron period 
27r/57c. Then in the term containing fcn performing integrations by parts in the w-integral, one finally obtains 
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^ — — oo 



kp? [k^a) A [^p (k,c.)] + |i- [Jf_, {k^a) - Jf+i (fc^a)] [ujP (k, . 
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The above expression is the main result of this paper For a simplicity we shall call the expression (|20] | as a Bessel-function 
representation of the stochastic energy loss rate. For many practical applications, however, it is important to represent the energy 
loss in an alternative but equivalent integral form, see Appendix lAl for details. In particular, the integral representation of the 
energy loss facilitates the limits of a heavy ion (with M ^ m, where m is the electron mass) and a weak magnetic field 
when is vanishingly small. Although the stochastic part of the energy loss, Eq. ( |20] |. vanishes in the limit M — oo for 
pedagogical purposes we treat this case first keeping the leading term only (proportional to M^^). In this case the projectile 
moves with rectilinear trajectory with an arbitrary angle with respect to the magnetic field. The limit fie of Sst is performed 
in AppendixlAl see Eq. ( IA3b . The second term of this expression contains a full derivative and vanishes after w-integration and 
the stochastic energy loss is then given by 

It should be noted that in Eq. ( |2TI ) the effect of a magnetic field is contained in the power spectrum P(k, w) of the fluctuations. 
It is seen that formally Eq. ( 1211 1 is equivalent to the stochastic energy loss derived in Refs. ||23| - |25| at vanishing magnetic field 
assuming that in Eq. (ISTT i P(k, w) is the power spectrum at _B = and neglecting the phase average J^^ dip /2tt. In the present 
context of the magnetized plasma this integration appears due to the average of the energy loss (l20t with respect to the cyclotron 
period iir/Vtc of the particle. Therefore for a conformity of our results with the case of unmagnetized plasma at i? the 
(y9-average is unnecessary and must be skipped. 

It is well known that un upper cutoff parameter fcmax = 1/ ''min (where rmin is the effective minimum impact parameter) must 
be introduced in Eqs. (fT4l i. (|20| | and (l2Tl i to avoid the logarithmic divergence at large k. This divergence corresponds to the 
incapability of the perturbation theory to treat close encounters between the projectile particle and the plasma electrons properly. 
For rmin we use the effective minimum impact parameter excluding hard Coulomb collisions with a scattering angle larger than 
7r/2. The resulting cutoff parameter /c^ax is well known for energy loss calculations (see, e.g., Refs. 10, [31 and references 
therein) and reads /cmax = ^{vq + ■i;^j^)/|Z|e^, where v^\^ — {T/mY/'^ is the thermal velocity of the electrons. 

Within linear response theory, the power spectrum P(k, w) of the fluctuations of the electrostatic fields at thermal equilibrium 
follows from the fluctuation-dissipation theorem and is completely determined by the longitudinal dielectric function e(k, of 
the medium lEllBil. 

AT -1 

P (k, uj) = ^Im^^ r. (22) 

{2TTfuj £(k,w) 

We consider magnetized electron plasma when the longitudinal dielectric function is given by (see, e.g., Ref. lEsIl ') 

e(k,t^) = l + -\^ Jl+ V K{z)[W{Pn)-l]\ , (23) 



n— — oo 



where z = k\a1, f3n — {lo — ncjc)/|fc|| |wth, A„(z) ~ e ^In{z), In is the modified Bessel function of the first kind, — Vth/ 



is the cyclotron radius of the electrons. Also W(C,) — g{C) + if{C) is the plasma dispersion function with the real (g(C)) and 
imaginary (/(C)) parts IH] 



In the next sections we consider some particular cases of the stochastic energy loss of the projectile particle, in particular at 
vanishing and infinitely strong magnetic fields. At vanishing magnetic field we denote the dielectric function as £o(fc, ^) which 
is isotropic and is given by the usual expression (see, e.g., Ref. |28]) 

eo{k,u;) = l + -^w(-^) . (25) 



To the dielectric function ( |25] | corresponds the isotropic power spectrum Po{k, u) according to Eq. {22\ . In this case the stochas- 
tic energy loss rate (1211 is simplified to ll23l - l25ll 

^"'"~M~ J Po{k,kvo)k^dk. (26) 

At infinitely strong magnetic field we denote the dielectric function as eoo(k,aj) which is determined by Eq. ( |25] l, where, 
however, the phase velocity w/fc in the argument of the plasma dispersion function is replaced by a;/|fc|| | (see, e.g., Ref. 
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The above expressions gives the mean energy (per unit time) absorbed (or emitted) by a propagating projectile particle from 
the heat reservoir. Physically, this arises from the energy absorption from thermal fluctuations. In the absence of the magnetic 
field it is to be noted here that because the spectral density Po(^,w) of the fluctuations of the electric fields is positive for 
positive frequencies by definition, according to Eq. (l26t the particle energy will grow due to interactions with the fluctuating 
electric fields. The present situation with a magnetized plasma may be quit different. It will be further discussed in Sec. Hill 



III. ELR WITHOUT AND WITH STRONG MAGNETIC FIELD 

To illustrate the problem of the stochastic energy loss and the interrelation of the treatments determined by Eqs. (Il4l and ( l20l ). 
we consider the cases without and with an infinitely strong magnetic field. In the absence of a magnetic field the interaction of 
the test particle with stochastic plasma electric fields yields a positive ELR (which corresponds to an energy gain by the particle) 
and the total ELR decreases. At strong magnetic field, however, the stochastic interaction of the projectile may result in positive 
or negative ELR depending on the velocity of the particle. 



A. ELRs at vanishing magnetic field 



For a vanishing magnetic field when ^ Ad (or ujc ^ <^p), the power spectrum of the fluctuation of the plasma electric 
fields is determined by Eq. (l22l l. where the longitudinal dielectric function is given by the field-free expression (l25i . At vanishing 
magnetic field the ELRs 5*0 and S'st have been evaluated previously in Refs. | f33l[34ll and ll23l - l25ll (see also references therein), 
respectively. In the leading order with respect to the cutoff parameter fcmax the ELRs Sq and S'st read 



So 



Ssi 




-/(A) 



S (u) udu 



^o^3^[/(A)lne-S(A)], 



(27) 
(28) 



where erf(z) is the error function, Qq = ArieZ'^e^ /mvth, £, — fcinaxAo = (A^ + l)/e\Z\, A = vo/vth is the projectile velocity 
in units of thermal velocity. Also 



E{u)^\f{u)\n [g' {u) + f {u)] + (u) 



arctan 



fin) 



(29) 



It is well known that at vanishing magnetic field the ELRs determined by Sq and S'st are always negative and positive, respec- 
tively. This features correspond to the energy loss and energy gain, respectively, by projectile particle. 

Consider briefly some limiting cases for the low- and high-velocity projectiles. Taking into account the behavior of the 
function ~ (l/2)^7r/2u (at u < 1) and E{u) ~ --njlu^ (at u > 1), at low- (A < 1) and high-velocities (A > 1) the 
ELRs become 



^0 
^0 



\' hi 



1 



Sst 



'211 In 



5st ^ Q^M^, 



1 



(30) 
(31) 



respectively. For derivation of the first relation in Eq. (1311 we have used the frequency sum rule for the dielectric function dZSl l. 
From Eqs. ( l30l l and dSTT l it is seen that for light projectile particle (for instance for an electron with M — m) the stochastic 
ELR may essentially contribute to the total energy loss at small velocities. Then the total rate So + S'st becomes positive which 
corresponds to an energy gain by projectile particle. At high-velocities, however, the ELR S'o essentially exceeds the stochastic 
energy loss and the total ELR remains negative. 



B. ELRs at infinitely strong magnetic field 

We now turn to the case when a projectile particle moves in a plasma with a strong magnetic field, which is, on one hand, 
sufficiently weak to allow a classical description (hujc < T), and, on the other hand, comparatively strong so that the cyclotron 
frequency of the plasma electrons exceeds the plasma frequency Wc ^ Wp (or fle ^ Ad)- This limits the values of the magnetic 
field, the temperature, and the plasma density. From these conditions we can obtain 3 x lO'^y^n^ < B < lO^T, where is 
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vJ10'cm/s] v^ii[10 cm/s] 

FIG. 1: (Color online) Electron energy loss rate (in units of 10'' MeV/s) in a plasma with rig = 10^" cm"'', T = 10 eV (e = 0.06) in terms of 
the velocity (in units of 10** cm/s) at vanishing (a) and strong (b) magnetic field. The energy loss rates 5*0 and Sst and the total rate So + Sat 
are shown as dashed, dotted and solid lines, respectively. 



measured in cm^'^, T in eV, and B in kG. The plasma electrons in the lowest order can respond to the excited waves only in the 
direction of B. Because of this assurnntion, the perpendicular cyclotron motion of the test and plasma particles is neglected. As 
shown by Rostoker and Rosenbluth ||5|] the plasma waves involved in this case are oblique plasma waves having the approximate 
dispersion relation cj — utpk^^/k. 

The dielectric function eoo(k, w) is obtained from Eq. dZSl l. where, however, the phase velocity Lu/k in the argument of the 
plasma dispersion function is replaced by w/|fc|||. The corresponding power spectrum of the fluctuations in the presence of 
strong magnetic field is connected to Soo (k, w) via relation ( |22] |. Neglecting the transversal cyclotron motion of the test particle 
and the electrons in Eqs. ( fT4l i and (|20] i in the lowest order one obtains (see, e.g., Ref. J^l^) 

^o--QoA[/(A)lne-S(A)], (32) 



>S'st — 



M 



where now A — VQ\^/vth- Substituting the power spectrum of the fluctuations with the dielectric function £00 (k, uj) in Eq. 
in the leading order with respect to the cutoff parameter ^ we arrive at 

777 (J 

^st = Qo^^[/(A)lnC-5(A)], (34) 

where S(A) is determined by Eq. ( |29] l. 

Consider briefly the limits of low- and high-velocities. At A ^ 1 from Eqs. (|32] | and (l34l i one obtains 

So^-Qoxf^X^(\n^-l-], (35) 



2 V 2 

^.t^Qo^./^fln^-iL (36) 



M V 2 V e\Z\ 2 
At high-velocities, A 1, Eqs. $32\ and ( l34b become 

^^^^ 

respectively. It is seen that unlike the field-free case where the energy loss rate Sst in Eqs. ( l30l l and (ISTT i is always positive, 
the ELR in the presence of a strong magnetic field is positive at small-velocities while it is negative at high-velocities. For a 
light particle (e.g. for an electron) the stochastic ELR even essentially exceeds the energy loss Sq at small velocity limit and the 
total ELR becomes positive. The results of the numerical evaluation of Eqs. ( l27b . ( l28T l and (|32] |. ( |34] | are shown in Fig. [T]for an 
electron test particle (AI = m) moving in a plasma with rig = 10^" cm""^ and T = 10 eV (e = 0.06). This figure shows the ELR 
(in units of 10^ MeV/s) as a function of the particle velocity (in units of 10^ cm/s) at vanishing (left panel) and at strong (right 
panel) magnetic field. It is evident that the effect of the electric field fluctuations on the test particle energy loss is significant at 
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low-velocities both for vanishing and strong magnetic field. At higher velocities the relative importance of the stochastic energy 
loss to the usual energy loss 5*0 decreases gradually. We also mention that similar calculation can be done for any test ion with 
M ^ m. However, the stochastic energy loss for an ion both at vanishing and infinitely strong magnetic field at least three 
orders of magnitude smaller than Sst shown in Fig.[T] 

We would like to close this section with the following remark. From Eqs. dSOl l and ( |36] | it is seen that the fluctuation energy loss 
rate does not vanish at low-velocities (as it occurs with the polarization energy loss rates d30t and d35t ) for vanishing and strong 
magnetic field and the corresponding friction coefficient diverges as vq 0. For arbitrary magnetic field similar conclusion 
can be made on the basis of a general expression (IA2t where the stochastic ELR is given in the integral representation. At the 
vanishing ion velocity (i.e. at vq 0) this expression yields 



Ost — 



M 



(38) 



It is seen that in low-velocity limit and for arbitrary B the rate Sst is finite and its longitudinal and transversal parts are deter- 
mined by static P{k, 0) and dynamic P{k, ftc) fluctuation spectra, respectively. At vanishing magnetic field (i.e. at 51c 0) 
the second term in Eq. (|38] | can be collected with the first one and the resulting ELR coincides with Eq. (|26] | in the limit vq 0. 
In this case the stochastic energy loss is only determined by the static fluctuation spectra P{k, 0). This feature of the stochastic 
ELR at vanishing velocity is closely resembled the Brownian motion. In fact at Mv^ /2 < T when the kinetic energy of the 
particle is comparable with the thermal energy of the plasma environment the particle motion is stochastic in nature. The energy 
change of the particle is determined by the correlation function C{t) — (vst(i) • Fst(i)), where Fst(i) is the stochastic force 
acting on the particle and Vgt (i ) is the stochastic velocity of the particle under the influence of Fgt [t). At low velocities Vgt {t) 
can be estimated employing the Langevin equation and the result at long times reads Vst(i) — {l/^M)Yst{t), where 7 is the 
friction coefficient. Thus, asymptotically C{t) ~ {l/^M)l^1^{t)) is a constant at thermal equilibrium which is determined 
by the thermodynamic properties of the medium. At vanishing magnetic field this constant is related to the static fluctuation 
spectrum 0) while in i? 7^ case the fluctuations give rise the cyclotron motion of the particle and the transversal part of 
C{t) is determined by the dynamic spectrum P(k, Vic), see the second term in Eq. 



IV. ELR AT LOW-VELOCITIES AND FOR INTERMEDIATE MAGNETIC FIELD 

In SecHni the effect of the electric-field fluctuations on the energy loss of the projectile particle has been investigated in two 
extreme regimes with vanishing and very strong magnetic fields. In this section we extend these results considering briefly some 
intermediate cases with finite magnetic field which cover the entire domain between field-free and strong magnetic field regimes. 
As will be shown below at ^ Ad the effect of the magnetic field on the energy loss process shows some additional features. 
In general, the analytical or numerical investigation of the energy loss rates given by Eqs. (fl4] i and (l20t with Eqs. (l22t and (l23T l 
for arbitrary B and the test particle velocity w is a formidable task and requires separate investigations. Previously the energy 
loss Sq has been investigated numerically and analytically for arbitrary B and v but mostly neglecting the gyration of the test 
particle, Refs. [2-16]. To gain more insight into the magnetic field effect on the stochastic energy loss process we consider here 
the limit of low-velocity projectile ion when electric field fluctuations are expected to be very important. In this case assuming a 
heavy ion {M ^ m) the imaginary part of the inverse dielectric function (spectral function) involved in Eq. (|22] | is represented 
in the form 



-1 Im[£(k, w)] Im[£(k,CL;)] 
"£(k:^ " |e(k,^)|^ ^ |£(k,0)r 



Im-7T— T = „ \ 2 2 ■ (39) 



Note that in the last part of Eq. (|39] l the dynamical screening has been replaced by the static one thus neglecting the collective 
plasma resonances involved in e(k, a;). For arbitrary velocities this simplified version of the spectral function and the corre- 
sponding ELR has been calculated in Refs. [2, 1 1] with vanishing and strong magnetic fields. It has been shown that at B — 
the simplified and the full spectral functions yield close ELRs at low-velocities. However, at t^o > wth when the collective 
effects play an important role these approaches deviate from each other On the other hand, the strong magnetic field reduces the 
transversal dynamics of the plasma electrons and the role of the collective effects becomes less important which results in the 
good agreement between two approaches in a whole velocity range. 

Consider the ELR Sq given by Eq. ( fT4] i assuming that the projectile ion moves along the magnetic field (i.e uo_l = 0) with 
I'd II = Vq <^ Vth- Inserting the approximate spectral function (|39) into Eq. ( fT4b and introducing new integration variables 
K = fc^Cg, fi = cos^ 9 (where 6 is the angle between k and B) we arrive at 

So^-.[^Qo\' f'df, g A„(^(l-;.))expf-;;^V (40) 
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where X = vq/ vth, rj = ag/Ao, and the other quantities have been introduced in Sees. lU and Hill 

The stochastic ELR at low-velocities is determined by Eq. dSST l where the first and the second terms are the contributions 
from longitudinal and transversal fluctuations, respectively. It should be emphasized that the relation ( |39] | is justified for the first 
term of Eq. ( |38] | where only the static spectrum P(k, 0) is involved. The second term contains the cyclotron frequency ftc of the 
projectile ion. Since ilc <C ojc the ion cannot resonate with plasma electron resonances involved in the dielectric function and in 
the fluctuation spectrum P(k, flc) the simplification given by Eq. ( l39b is again justified. 

Consider now Eq. ( [38] l with Eqs. (l22T i and (|39] l. Using the integration variables introduced above we obtain 



{n — a)^ 



X { /^exp ( + (1 - At)exp 



2KfJ. 



where a = \Z\m/M ^ 1. Let us recall that at low-velocities the ELR, Sst, in the leading order does not depend on the ion 
velocity while 5*0 ^ (see Eq. (|40]|). It is seen that the term of Eq. (1411 1 with n = (no cyclotron motion) exhibits a singular 
behavior in the limit of a ^ 1 where the ^ integration diverges logarithmically for small ji (or small fc||). We must therefore 
keep a finite in that integration to avoid such a divergence. This anomalous contribution that arises from the term of Eq. (|411 
with n = is denoted as S'^^. The other terms of Eq. (l4ll with n 7^ do not exhibit singular behavior at a ^ 1 and their 
contribution to the stochastic energy loss is denoted as S'gj {Sst ~ -Sgt + "^^st)- Thus in the limit of a heavy projectile ion with 
a ^ 1 from Eq. W\\ in the leading order we obtain 



2 c2 

IT m 



« KdK r , , , /, 2k 
M Jo (n + 'q^y [ V a 



SL-d^i,Qol Ao(.) ln9-7) (42) 



1 



+ / [Ao{K{l-^i))-Ao{K)]^ 
Jo A* 

- r 7^ - £ (1 - (-f-) ' ^''^ 

V2M Jo (K + Tj^yJo M V 2KAi/ 

where 7 = 0.5772 ... is the Euler's constant. The terms neglected in deriving Eqs. ( |42] | and (|43] | are of the order of O(a^) and 
0{a^ In a). Also we note that the anomalous ELR S'^^ given by Eq. (l42l i is valid at a < 77 (or 17^ < Wp) which is easily fulfilled 
for heavy ions. In the opposite case with a > rj (but a 1) one should consider more general Eq. (1411 as a starting point for 
investigation of the anomalous ELR. 

The ELR ( |42] | can be further simplified assuming that 77 > 1/^ (or Oe > rmin)- This condition allows one to extend the upper 
limit of the integration in Eq. (l42T i to infinity which yields 



TT m 
8 m' 



(44) 



where 



Go (x) = / °^ \, = [(1 + x) Ko (x) - xK, {x)] , (45) 
Jo (k + x) 

Gi(x)= / In (k + x) ^ ^ ■ (46) 

Jo {k + x) 

Here (x) and Ki (x) are the modified Bessel functions of the second kind. The integral involved in the function Go (x) has 
been evaluated by Ichimaru (see Ref. ||28|] for details). 

Consider briefly the connection of Eqs. (l4Qt-(l44t to the ELRs obtained in the previous section for two limiting cases. It is 
seen from Eqs. (|40] | and (ITTT i that at very strong magnetic field only the terms with n = contribute to the ELRs. Then the 
ELRs ( |40l l and iAli in the limit of 77 — > after k integration coincide with Eqs. ( [35] l and ( |36] l. respectively. We merely note that 
in this limit the ELR 5"^ vanishes and the stochastic energy loss in the case of very strong magnetic field is determined by the 
anomalous term S'gj. 

Next we treat the limit of the vanishing magnetic field (77 00). In this case all cyclotron harmonics in Eqs. (|40] | and dTIT l 
contribute to the ELRs. It is useful to transform the 77, summations in these expressions into an integral form where the limit 
77 — > (X) can be easily done. Since this procedure is similar to one considered in Appendix lAl we refer a reader to Refs. Ij2j 1I2I 
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B/B 



FIG. 2: (Color online) Proton energy loss rate (in units of 10^ MeV/s) in a plasma with rie — 10^" cm~^, T = 10 eV (e = 0.06) vs scaled 
magnetic field B/Bs and at low-velocities (A = 0.1). The energy loss rates 5*0, S'^t and Sst and the total rate So + Sat are shown as dashed, 
dot-dashed, dotted and solid lines, respectively. See the text for explanations. 

for more details. As a result we obtain that at vanishing magnetic field the n summation in Eq. (l40l l is replaced by ^/JI. The 
remaining integrals can be easily performed. After evaluating these integrals one can see that the ELR Sq at 77 — > oo completely 
agrees with the first relation in Eq. d30i . We turn now to the anomalous term given by Eq. (l44l l. Using the asymptotic behavior of 
the modified Bessel functions one cane see that the functions Go(?7^) and Gi{rf') at ?/ — )• cx) behave as ^ I/77 and ^ (I/??) In 77, 
respectively, and thus vanish in the field-free case. On the other hand the stochastic ELR S'^J. is evaluated in a same manner as 
explained above and the summation in Eq. (|43] | in a limit of vanishing magnetic field yields ^JJH/'i. Substituting this formula 
into Eq. ( |43] l and performing ji and k, integrations we finally arrive at the second relation of Eq. (|30l l. The anomalous term 
Eqs. (l42b and ( l44l l therefore represent a new effect arising from the presence of the magnetic field. The physical origin of such 
an anomalous ELR is the same as those obtained in Ref. |28] and Refs. |2, 10] for the transport coefficients and the stopping 
power of an ion, respectively. It may be traced to the spiral motion of the electrons along the magnetic field lines. These electrons 
naturally tend to couple strongly with long-wavelength fluctuations (i.e., small fc||) along the magnetic field. In addition, when 
such fluctuations are characterized by slow variation in time (i.e., small u) — ilc), the contact time or the rate of energy exchange 
between the electrons and the fluctuations will be further enhanced. In a plasma, such low-frequency fluctuations are provided 
by the slow projectile ion. The above coupling can therefore be an efficient mechanism of energy exchange between the electrons 
and the projectile ion. In the limit of fie (or M — ^ cjo), the contact time thus becomes infinite and the ELR diverges. 

In the presence of strong magnetic field one could expect that for heavy ions the anomalous ELR S'^^ mainly contributes to 
the total stochastic ELR 5st. To demonstrate this feature the ELRs are plotted in Fig. |2]vs scaled magnetic field B/Bg for 
Vf) = O.lwth — 1-3 X 10'' cm/s, — 10^" cm^'^, T = 10 eV, where Bs = mavp/e ~ 3.2 x lO'' kG. The energy loss rates 
So, /Sgt and ^st and the total rate Sq + Sst are found for a proton projectile and are shown as dashed, dot-dashed, dotted and 
solid lines, respectively. In Fig.|2]the limiting cases of the vanishing and very strong magnetic fields considered in Sec. Hill are 
reached at B < Bs and B > lO^Bg, respectively. It is seen that the magnetic field may essentially increase the ELRs in plasma 
compared to the field-free (Eq. (|30] l) and strong field (Eqs. (|35] l and ( |36] l) regimes and this occurs at B > lOBg. In the latter 
case the stochastic ELR is mainly determined by the anomalous ELR S!.^. Furthermore, at low ion velocities the fluctuations in 
the electric field essentially counter the effect of the polarization field resulting in the reduced ELR (solid line in Fig.|2]i. 

As discussed above, Eqs. (|40] i and (HTb and hence the results shown in Fig. |2] account for only static screening neglecting 
dynamic polarization effects. As shown in Refs. ^ [TTt] in a vanishing or very strong magnetic field these effects are not 
important at low-velocities, and the situation with finite B requires further investigations, in particular for electron projectiles. 
Intuitively it is clear that a gyrating test electron may resonate with a electron plasma fluctuations with a cyclotron frequency Wc, 
and the approximation ( [39] l is not as obvious as for the case of a heavy ion. Also it should be emphasized that the validity of the 
regime of a classically strong magnetic field requires the condition B < Be = {mc/eK)T. Thus the results shown in Fig.|2]are 
valid up to Bc/Bs = T/fiuip ~ 27. Clearly the realization of the regime of a strong magnetic field requires high temperatures 
and low densities and the enhancement of the ELRs at i? > 10i?s may not be accessible under certain conditions. 

V. SUMMARY 

In this paper within linear response theory we have investigated the energy loss rate (ELR) of a nonrelativistic projectile 
particle moving in a magnetized classical plasma due to the electric field fluctuations. In the course of this study we have 
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derived general expression for the stochastic energy loss rate. As in the field-free case the later is completely determined by the 
fluctuation power spectrum P(k, u). Assuming an equilibrated and homogeneous magnetized plasma when the power spectrum 
of the fluctuations is derived from the fluctuation-dissipation theorem and is expressed via longitudinal dielectric function of the 
medium we have considered two somewhat distinct cases with vanishing and with extremely strong magnetic field. At vanishing 
magnetic field the electric field fluctuations lead to an energy gain of the charged particle for all velocities. Physically, this is 
arisen due to the energy absorption by projectile particle from thermal fluctuations. It has been shown that in the presence of 
strong magnetic field this effect occurs only at low-velocities. At high-velocities the test particle systematically loses its energy 
due to the interaction with a stochastic electric field. 

In the low-velocity limit but for arbitrary magnetic field, we have found an enhanced stochastic ELR of an ion compared to 
the limiting cases with vanishing and very strong magnetic fields, mainly due to the strong coupling between the spiral motion 
of the electrons and the long-wavelength, low-frequency fluctuations excited by the projectile. The nature of this anomalous 
ELR is only conditioned by the external magnetic field. 

The analysis presented above can in principle be extended to the case of a nonequilibrium magnetized plasma, provided the 
power spectrum of the electrostatic fluctuations are known. In principle the power spectrum can be derived within linear response 
theory lK4ll25ll . Fluctuations are much stronger in nonequilibrium situations than in systems in thermal equilibrium [24, 25']. 
The effect of the electric field fluctuations on the particle energy loss is therefore expected to be stronger in a nonequilibrium 
plasma. We intend to address this issue in our forthcoming investigations. 
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Appendix A: Integral representation of tlie stochastic energy loss 

In Sec. [n] we have derived the Bessel-function representation of the stochastic energy loss of a gyrating test particle. For 
some applications an integral representation of Sst is desirable. For deriving the integral form of the energy loss we rewrite the 
delta- function S{uj — Q (k)) in Eq. (|20] i using an integral representation of this function. Then Eq. (|20] i becomes 



du die^('^-'=ii"oii)t e-''"<=* Jf (fc^a) (Al) 

^ J — OO J — CO ]_ 



1— — OQ 



Then using the summation formula J2'iZ-oo^^''*'^ii'^) ~ «/o(2zsin|) 1.3 1.1 the stochastic energy loss may be alternatively 
represented in the form 



Sst = I I I -^0 ( 2fc_Lasin^ ) dt (A2) 



d sin (Vict) 

kf\-g^ [wP (k, uj)] cos [{uj - fc||t;o||) t] + [ujP (k,a;)] — ^^-^ sin [{u - k\\v^\^ t] 

Using this integral form it is now easy to perform the limit of a heavy ion (with M — > oo) when fie 0. In this limit the zero 
order Bessel function becomes jQ{k±VQ±t) which we represent in the integral form as an average of the exponential function 
^±ik±-vo±t ^jjjj j-gspgct to the azimuthal angle of the vector vo±. Thus performing ^-integration the energy loss (IA2t in the 
limit fie — >^ becomes 

s^.ii^riif^r^ (A3) 

— oo 

2 



M /n 2n 



<5 (c^ - k • vo) ^ [o.F (k, u)]^'j^-^[6{uj~k- vo) ojP (k, u)] 
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The second term in Eq. (IA3l l contains full derivative and vanishes after w-integration. Then the averaged stochastic energy loss 
is solely given by the first term, see Eq. (ISTT l. 
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